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If dark matter particles are not perfectly stable, their decay products might be seen in the cosmic- 
ray fluxes. A natural candidate for decaying dark matter is the gravitino in R-parity violating 
scenarios. In the relevant GeV-TeV energy range, the Fermi Large Area Telescope (LAT) is now 
measuring cosmic gamma-ray fluxes with an unprecedented precision. We use the Fermi LAT 
gamma-ray data to search for signatures from gravitino dark matter particles, concentrating on 
gamma-ray lines and galaxy cluster observations. Implications of our results for the decay length of 
the next-to-lightest superparticle, which could be seen at the LHC in the near future, are discussed. 
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I. INTRODUCTION 

A theoretically well motivated example for decay- 
ing dark matter is the gravitino 1P3/2, which appears 
in locally supersymmetric extensions of the Standard 
Model. In scenarios where -R-parity is mildly vio- 
lated and the gravitino is the lightest superparticle 
(LSP), thermal leptogenesis, gravitino dark matter 
and primordial nucleosynthesis are naturally consis- 
tent Within this framework, the gravitino would 
decay with cosmological lifetimes 0], making its de- 
cay products potentially observable in the cosmic-ray 
fluxes • For gravitino masses m 3 / 2 < 200 GeV, 
the most prominent feature in the decay spectrum is 
an intense gamma-ray line, produced by two-body de- 
cay into neutrinos and photons, 1P3/2 —> iv 0|; this 
line can be searched for in the gamma-ray fluxes ob- 
served at high latitudes. For larger gravitino masses 
m 3/2 <; 200 GeV, the branching ratio into gamma- 
ray lines is suppressed, and instead the decay modes 
^3/2 - ► W ± £ T and ^3/2 — ► produce a gamma- 
ray flux with a broad continuous energy spectrum; 
this flux could potentially show up in observations of 
galaxy clusters or the extragalactic gamma-ray back- 
ground (see e.g. Ref. [3, H, 0)- 

Here, we briefly summarize our searches for gamma- 
ray signals from gravitino dark matter in the data 
of the Fermi Large Area Telescope (LAT) [lf3 |. 
Firstly [Til ]. we extend the gamma-ray line analysis 
presented in Ref. [T^j to a larger energy range of 1- 
300 GeV, searching for significant line signals that 
might come from dark matter decay (or annihilation) 
in the Galactic dark matter halo. Secondly [l3| . we 
analyze the gamma-ray flux from eight galaxy clus- 
ters, targets which are more sensitive to continuous 
spectra. Extending previous analysis (l^.[l5|. we treat 
the decaying dark matter signal as extended source 
and analyze the different target clusters individually 



as well as in a combined likelihood approach. We 
present constraints on the dark matter lifetime as well 
as on the annihilation cross section. We then apply 
our findings to the scenario of decaying gravitino dark 
matter and comment on implications for the possible 
observation of long-lived superparticles at the LHC. 

The remaining sections are organized as follows: In 
the second section, we introduce briefly the gravitino 
dark matter scenario, in section three we summarize 
our gamma-ray line and galaxy cluster analysis, and 
in the fourth section we present the resulting limits 
on the gravitino lifetime and decay width of the next- 
to- lightest superparticle (NLSP). 



II. GRAVITINO DARK MATTER 

Among the different scenarios that were proposed 
to reconcile thermal leptogenesis and gravitino dark 
matter with the standard BBN scenario [l6l - [l9j . a 
mild violation of .R-parity that induces a rapid decay 
of the NLSP before the onset of the BBN is maybe 
the most interesting from the perspective of indirect 
dark matter searches 0. If imparity is violated, the 
gravitino dark matter particle becomes unstable and 
subject to decay, which opens the possibility to look 
for its decay products in the cosmic-ray fluxes. 

Here, we consider the supersymmetric standard 
model with explicit bilinear .R-parity violation as de- 
scribed in Ref. [201 ] . Trading the mass mixing param- 
eters for -R-parity breaking Yukawa couplings as pro- 
posed in Ref. [2fJ] , the gravitino decay is a function of 
a single dimensionless parameter £, which also enters 
the decay width of the NLSP. As boundary condi- 
tions for the supersymmetry breaking parameters of 
the MSSM at the grand unification (GUT) scale, we 
consider equal scalar and gaugino masses mo = m 1( / 2 , 
a zero trilinear scalar coupling ao — 0, and tan /3 = 10. 
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FIG. 1: Two-body decay branching ratios of the gravitino, 
from Ref. [Hj]. 



In this case, the bino-like neutralino Xi is the NLSP; 
the universal gaugino mass m 1 / 2 remains as the only 
independent variable, and the gaugino masses 
satisfy the following relations at the electroweak scale: 
M 3 /Mi ~ 5.9 and M-JMa ~ 1.9. Electroweak pre- 
cision tests (EWPT) |2CH22|], as well as the possible 
overproduction of gravitinos in presence of the high 
reheating temperatures required by standard thermal 
leptogenesis [23] , yield further bounds on the gravitino 
mass like TO3/2 > 30 GeV, and on the NLSP neutralino 
mass like 100 GeV < m x o < 690 GeV (for details see 
Hef. [LH). 

The gravitino inverse deca y rate into pho- 
ton/neutrino pairs is given by [2I |20|| 

x _ 32\/2 GpMp MfM| 

V-3/2->7" ~ m 3 ( Ma _ Ml )2 ' W 

where a is the electromagnetic fine structure constant, 
M P = 2.4 x 10 18 GeV the reduced Planck mass, and 
G F = 1.16 x 10~ 5 GeV~ 2 is the Fermi constant. The 
branching ratios into other channels are presented in 
FiglU 



III. FERMI LAT LIMITS 
A. Gamma- Ray Lines 

Our line search is based on the measurements of the 
cosmic gamma-ray flux performed by the Large Area 
Telescope (LAT). The gamma-ray events that enter 
our analysis are selected from the 'DataClean' event 
class measured between 4 Aug 2008 and 17 Nov 2010. 
We consider energies between 1 GeV and 300 GeV, 
and apply the zenith angle criterion 9 < 105° in or- 
der to avoid contamination by the Earth's Albedo. 
The expected shape of the measured gamma-ray line 
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FIG. 2: Lower bounds on the dark matter inverse de- 
cay width into monochromatic photons and neutrinos as 
a function of the dark matter mass m^,, cp. Ref. The 
gray-solid line shows the 95% C.L. limits as function of 
the gamma-ray line energy, the black dots show the weak- 
est limits obtained in certain adopted energy bands. The 
previous Fermi LAT limits from Ref. are shown for 
comparison. 

spectrum is inferred from the Fermi LAT instru- 
ment response function. In our analysis, we consider 
only Galactic contributions to the dark matter sig- 
nal, and take the Navarro-Frenk- White (NFW) pro- 
file [ill, [24j as a reference for the dark matter distri- 
bution (Einasto or isothermal profiles would lead to 
very similar results). All profiles are normalized to 
Pdm = 0.4 GeV cm -3 at Sun's position. For decaying 
dark matter signals we choose to consider the whole 
sky excluding only the Galactic disk at |6| < 10° with 
its large foregrounds, since this large region features 
the best signal-to-noise ratio. 

The profile likelihood method [25| is used to calcu- 
late the significance of a potential gamma-ray line con- 
tribution to the observed gamma-ray flux. The data 
are modeled by a simple power law plus a line signal at 
fixed energy E-y . Since the power law is only locally a 
good approximation to the background fluxes, we use 
a small sliding energy window in the fitting procedure. 
The size of this energy window varies between ±2a4 8 jf 
at low gamma-ray line energies E 1 , and roughly 07 7 
to 3E 7 at high gamma-ray line energies. Lifetime up- 
per limits at the 95% C.L. are derived by increasing 
the line signal and refitting the remaining parameters 
until the -21og(likelihood) of the fit increases by 4 from 
its best-fit value. 

No gamma-ray lines with 5<r significance were found 
in our analysis; the corresponding limits on dark mat- 
ter decay into monochromatic photons are shown in 

Fig. La 



B. Galaxy Cluster Observations 

The eight galaxy clusters that we consider in this 
work are Fornax, Coma, A1367, A1060, AWM7, S636, 
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FIG. 3: Lower limits on dark matter lifetime for decay 
into bb final states, as function of the dark matter mass, 
cp. Ref. [13] . Solid lines show individual cluster limits, 
the dashed line the limits from the combined likelihood 
analysis. The dotted line shows for comparison the limit 
that can be derived from the EGBG. 



NGC4636 and NGC5813. They are selected from the 
extended HIFLUGCS X-ray catalog [H S3 in or- 
der to yield large signals from dark matter decay. 
The gamma-ray events entering our analysis are taken 
from the P7S0URCE_V6 event class of the Fermi LAT 
data measured between 4 Aug 2008 and 21 Jul 2011. 
From all events recorded by the Fermi LAT, we select 
those with energies between 400 MeV and 100 GeV 
and apply the zenith angle criterion 9 < 100° in order 
to avoid contamination by the Earth's Albedo. For 
each galaxy cluster, we consider photons events in a 
10° x 10° squared region centered on the cluster posi- 
tion. These events are binned into a cube of 0.1° x 0.1° 
pixels with 24 logarithmic energy bins. 

We assume that the smooth component of the dark 
matter halo follows a Navarro-Frenk- White (NFW) 
profile [HI, [24j, where the scale radius r s and the 
density normalization p s have to be determined from 
observations. We adopt the observationally obtained 
concentration-mass relation from Ref. [28[ and use the 
cluster masses derived from ROSAT PSPC X-ray ob- 
servations in the extended HIFLUGCS catalog [27j to 
calculate the signal surface densities. Gamma rays 
from inverse Compton scattering between the CMB 
and the electrons and positrons that are produced in 
the dark matter decay are fully taken into account. 

As above, we use the profile likelihood method to fit 
the data with background and signal fluxes |25|. For 
the diffuse background fluxes we take the isotropic 
emission and the galactic foreground model templates 
currently advocated by the Fermi LAT collabora- 
tion for point source analysis (iso_p7v6source and 
gal_2yearp7v6_v0). On top of the diffuse templates, 
we add the point sources from the second Fermi LAT 
catalog 2FGL [29| within a radius of 12° around the 
cluster centers, as well as the extended dark matter 
signal. We account for uncertainties of the cluster 




FIG. 4: Lower limits on the gravitino lifetime, cp. Ref. [131 ]. 
The dot-dashed line shows the gamma-ray line limits, the 
dashed line the limits resulting from the combined cluster 
analysis (using the branching ratios shown in Fig.[T]), and 
the dotted line the EGBG limits. 



masses as determined by X-ray observations as a sys- 
tematic errors in our analysis. To this end, we approx- 
imate the posterior probability for the cluster masses 
by log-normal distributions. The resulting signal un- 
certainties are as large as a factor two in some cases. 
Finally, to combine the statistical power of the dif- 
ferent target regions and to reduce the impact of the 
cluster mass uncertainties, we performed a combined 
likelihood analysis of all eight clusters simultaneously. 
In this case, the combined likelihood function is de- 
fined as the product of the likelihood functions for the 
individual clusters. The only parameter that is bound 
to be identical for all targets is the dark matter life- 
time. 

No significant emission from the target clusters was 
found. For the case of decay into bb, our resulting 
limits are shown in Fig. [3] There, we show the limits 
that we obtain from the clusters individually (solid 
lines), as well as the limit from the combined anal- 
ysis (dashed line). The dotted line shows limits de- 
rived from the extra-galactic gamma-ray background 
(EGBG) as measured by Fermi LAT in Ref. [30| (see 
Ref. [Rf for details). 



IV. DISCUSSION 

In Fig. [4] we finally summarize the limits on the 
gravitino lifetime that we obtain from our gamma-ray 
line searches (dotdashed) and the galaxy cluster anal- 
ysis (dashed). As one can see from this figure, while 
the search for gamma-ray lines is efficient for grav- 
itino masses m 3 / 2 < 200 GeV, constraints from galaxy 
clusters observations dominates for m 3 / 2 > 200 GeV. 
However, at high gravitino masses even stronger limits 
come from the measured EGBG (dotted). 

A neutralino NLSP heavier than 100 GeV domi- 
nantly decays into IU ± ^ =F and Z°i>. The correspond- 
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Neutralino NLSP: Minimal allowed decay length [m] 
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FIG. 5: Contour plot of lower bounds on the neutralino 
NLSP decay length, coming from both gamma-ray line and 
cluster constraints on the gravitino lifetime, as function of 
the neutralino and gravitino masses, m x o and m 3 / 2 re- 
spectively, cp. Ref. [l3l ]. The lower gray region is excluded 
by electroweak precision tests (EWPT). For thermal lep- 
togenesis, overproduction (O.P.) of gravitinos excludes at 
minimum the left green region, a limit which strengthens 
to the black-dashed line when assuming the discussed uni- 
versal boundary conditions. 



ing decay width is directly proportional to £ squared, 
which also enters the gravitino decay width Eq. (Q]). 
As a consequence, the two quantities can be directly 
related [2fj| . Using the above gaugino mass rela- 
tion, lower bounds on the neutralino decay length 
CT-oo can be derived. Our results are summarized in 
Fig. [5J For the parameter space allowed by EWPT 
and overproduction bounds, we obtain minimal decay 
lengths 0(100 m — 100km); decay lengths as small 
as ct^o ~ 60 m are allowed for m 3 / 2 — 30 GeV at 
m x o ~ 140 GeV. Most interestingly, even these large 
decay lengths are in the range of detectability of the 
LHC [5|, since processes like x? ~> Z Q v — > /i + [jTv 
with displaced vertices inside the collider _would lead 
to essentially background free signatures 3l!]. 
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